Fluid-structure interaction simulation of prosthetic aortic valves : comparison between immersed boundary and arbitrary Lagrangian-Eulerian techniques for the mesh representation

In recent years the role of FSI (fluid-structure interaction) simulations in the analysis of the fluid-mechanics of heart valves is becoming more and more important, being able to capture the interaction between the blood and both the surrounding biological tissues and the valve itself. When setting up an FSI simulation, several choices have to be made to select the most suitable approach for the case of interest: in particular, to simulate flexible leaflet cardiac valves, the type of discretization of the fluid domain is crucial, which can be described with an ALE (Arbitrary Lagrangian-Eulerian) or an Eulerian formulation. The majority of the reported 3D heart valve FSI simulations are performed with the Eulerian formulation, allowing for large deformations of the domains without compromising the quality of the fluid grid. Nevertheless, it is known that the ALE-FSI approach guarantees more accurate results at the interface between the solid and the fluid. The goal of this paper is to describe the same aortic valve model in the two cases, comparing the performances of an ALE-based FSI solution and an Eulerian-based FSI approach. After a first simplified 2D case, the aortic geometry was considered in a full 3D set-up. The model was kept as similar as possible in the two settings, to better compare the simulations' outcomes. Although for the 2D case the differences were unsubstantial, in our experience the performance of a full 3D ALE-FSI simulation was significantly limited by the technical problems and requirements inherent to the ALE formulation, mainly related to the mesh motion and deformation of the fluid domain. As a secondary outcome of this work, it is important to point out that the choice of the solver also influenced the reliability of the final results

Immersed boundary-finite element model of fluid-structure interaction in the aortic root

It has long been recognized that aortic root elasticity helps to ensure efficient aortic valve closure, but our understanding of the functional importance of the elasticity and geometry of the aortic root continues to evolve as increasingly detailed in vivo imaging data become available. Herein, we describe fluid-structure interaction models of the aortic root, including the aortic valve leaflets, the sinuses of Valsalva, the aortic annulus, and the sinotubular junction, that employ a version of Peskin's immersed boundary (IB) method with a finite element (FE) description of the structural elasticity. We develop both an idealized model of the root with three-fold symmetry of the aortic sinuses and valve leaflets, and a more realistic model that accounts for the differences in the sizes of the left, right, and noncoronary sinuses and corresponding valve cusps. As in earlier work, we use fiber-based models of the valve leaflets, but this study extends earlier IB models of the aortic root by employing incompressible hyperelastic models of the mechanics of the sinuses and ascending aorta using a constitutive law fit to experimental data from human aortic root tissue. In vivo pressure loading is accounted for by a backwards displacement method that determines the unloaded configurations of the root models. Our models yield realistic cardiac output at physiological pressures, with low transvalvular pressure differences during forward flow, minimal regurgitation during valve closure, and realistic pressure loads when the valve is closed during diastole. Further, results from high-resolution computations demonstrate that IB models of the aortic valve are able to produce essentially grid-converged dynamics at practical grid spacings for the high-Reynolds number flows of the aortic root

Numerical Investigation of the Performance of Three Hinge Designs of Bileaflet Mechanical Heart Valves

Thromboembolic complications (TECs) of bileaflet mechanical heart valves (BMHVs) are believed to be due to the nonphysiologic mechanical stresses imposed on blood elements by the hinge flows. Relating hinge flow features to design features is, therefore, essential to ultimately design BMHVs with lower TEC rates. This study aims at simulating the pulsatile three-dimensional hinge flows of three BMHVs and estimating the TEC potential associated with each hinge design. Hinge geometries are constructed from micro-computed tomography scans of BMHVs. Simulations are conducted using a Cartesian sharp-interface immersed-boundary methodology combined with a second-order accurate fractional-step method. Leaflet motion and flow boundary conditions are extracted from fluid–structure-interaction simulations of BMHV bulk flow. The numerical results are analyzed using a particle-tracking approach coupled with existing blood damage models. The gap width and, more importantly, the shape of the recess and leaflet are found to impact the flow distribution and TEC potential. Smooth, streamlined surfaces appear to be more favorable than sharp corners or sudden shape transitions. The developed framework will enable pragmatic and cost-efficient preclinical evaluation of BMHV prototypes prior to valve manufacturing. Application to a wide range of hinges with varying design parameters will eventually help in determining the optimal hinge design

Simulation of the Three-Dimensional Hinge Flow Fields of a Bileaflet Mechanical Heart Valve Under Aortic Conditions

Thromboembolic complications of bileaflet mechanical heart valves (BMHV) are believed to be due to detrimental stresses imposed on blood elements by the hinge flows. Characterization of these flows is thus crucial to identify the underlying causes for complications. In this study, we conduct three-dimensional pulsatile flow simulations through the hinge of a BMHV under aortic conditions. Hinge and leaflet geometries are reconstructed from the Micro-Computed Tomography scans of a BMHV. Simulations are conducted using a Cartesian sharp-interface immersed-boundary methodology combined with a second-order accurate fractional-step method. Physiologic flow boundary conditions and leaflet motion are extracted from the Fluid–Structure Interaction simulations of the bulk of the flow through a BMHV. Calculations reveal the presence, throughout the cardiac cycle, of flow patterns known to be detrimental to blood elements. Flow fields are characterized by: (1) complex systolic flows, with rotating structures and slow reverse flow pattern, and (2) two strong diastolic leakage jets accompanied by fast reverse flow at the hinge bottom. Elevated shear stresses, up to 1920 dyn/cm2 during systole and 6115 dyn/cm2 during diastole, are reported. This study underscores the need to conduct three-dimensional simulations throughout the cardiac cycle to fully characterize the complexity and thromboembolic potential of the hinge flows

The fish tail motion forms an attached leading edge vortex

We propose an approach to falsification of oscillation properties of parametric biological models, based on the recently developed techniques for testing continuous and hybrid systems. In this approach, an oscillation property can be specified using a hybrid automaton, which is then used to guide the exploration in the state and input spaces to search for the behaviors that do not satisfy the property. We illustrate the approach with the Laub-Loomis model for spontaneous oscillations during the aggregation stage of Dictyostelium.Comment: In Proceedings HSB 2013, arXiv:1308.572

Hydrodynamics of the bluegill sunfish C-start escape response: three-dimensional simulations and comparison with experimental data

In this work we study the hydrodynamics of a bluegill sunfish performing a C-start maneuver in unprecedented detail using 3-D numerical simulations guided by previous laboratory experiments with live fish. The 3-D fish body geometry and kinematics are reconstructed from the experiments using high-speed video and prescribed as input to the numerical simulation. The calculated instantaneous flow fields at various stages of the C-start maneuver are compared with the two-dimensional particle image velocimetry measurements, and are shown to capture essentially all flow features observed in the measurements with good quantitative accuracy; the simulations reveal the experimentally observed three primary jet flow patterns whose momentum time series are in very good agreement with the measured flow field. The simulations elucidate for the first time the complex 3-D structure of the wake during C-starts, revealing an intricate vortical structure consisting of multiple connected vortex loops at the end of the C-start. We also find that the force calculated based on the 3-D flow field has higher magnitudes than that implied by the jet momentum on the midplane, and it exhibits large and rapid fluctuations during the two stages of the C-start. These fluctuations are physical and are related to the change in the direction of the acceleration of the fish body, which changes the location of the high and low pressure pockets around the fish.Organismic and Evolutionary Biolog